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Abstract 



We consider the moment space Mn corresponding to p x p real or complex matrix 
measures defined on the interval [0, 1]. The asymptotic properties of the first k components 
of a uniformly distributed vector {Si^n, ■ ■ ■ , Sn,n)* ^ U{M.n) are studied if n ^ oo. In 
particular, it is shown that an appropriately centered and standardized version of the vector 
{Si^n-, ■ ■ ■ 1 Sk,n)* converges weakly to a vector of k independent p x p Gaussian ensembles. 
For the proof of our results we use some new relations between ordinary moments and 
canonical moments of matrix measures which are of own interest. In particular, it is shown 
that the first k canonical moments corresponding to the uniform distribution on the real or 
complex moment space Ain are independent multivariate Beta distributed random variables 
and that each of these random variables converge in distribution (if the parameters converge 
to infinity) to the Gaussian orthogonal ensemble or to the Gaussian unitary ensemble, 
respectively. 

Keyword and Phrases: Gaussian ensemble, matrix measures, canonical moments, multivariate 
Beta distribution, Jacobi ensemble, random matrix 



of signed real (complex) measures /Xjj, such that for each Borel set A C [0, 1] the matrix 
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1 Introduction 




matrix measure /i on the interval 
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/i(y4) = {j2ij{A))^ is symmetric (hermitian) and nonnegative definite. Additionally, we re- 
quire the matrix measure to be normalized, that is /x([0, 1]) = Ip, where Ip denotes the p x p 
identity matrix. In recent years considerable interest has been shown in generalizing many of 
the results on classical moment theory, orthogonal polynomials, quadrature for mulas e tc. to 
the case of matrix measures. Am ong many others w e refer to the early pape r of iKreiru (119491) 



fll990r). ISinap and van Assche (|l994h. iDuran and van Asschel (Il99il 



and IPuran and Lopez- Rodrigue Jl996l . Il997h . ICriinbauml toO^ 



Duran 



and t o the more recent work s of Geronimol (1982[). Aptekarev and Nikishin _( 1983h 



19951. 



Rodman 



19961. 



I999I ) 



Griinbaum et al. 



([20051) and 



Damanik et al.l (120081 ) among many others. 

The aim of the present paper is to explore the relations between moments of matrix m easures 
and G aussian ensembles, an important distribution in the area of random matrices [see iMehta 
( I2OO4J )]. Both fields have been investigated rather independently and in this paper we demon- 
strate that there exists a deep connection between random moments and Gaussian ensembles, if 
the "dimension" of the moment space converges to infinity. To be precise, consider the real case 
and recall that the moments of a real matrix measure /z on the interval [0, 1] are defined by 



Sk= [ x^d^{x) e 5p(R); A; = 0, 1, 2, 
Jo 



and the nth moment space is given by 



(1.2) MnW = USl,...,Sn) 



x'dfiix), j = l,...,n} C (SpiR)) 



where Sp (M) denotes the set of all real symmetri c pxp matrices. In th e scala r case p = 1 this space 
has be e n investiga t ed by n umerous authors [se e 
Jl966h . ISkibinskvl Jl967h. lOette and Studded ( 



Karlin and Shapeleyl (119531 ) . iKarlin and Studden 



19971 )] and some of these r esults have been gen- 



eralized to the matrix case [see IChen and Li[(ll999h . IPette and Studded (l2002h among others]. 
In order to underst and the geometric properties of the moment space A^„(R) in the case p = 1, 
Chang et al.l (119931 ) proposed to consider a uniform distribution on A^„(R) and studied the 
asymptotic properties of random moment vectors. In particular, these authors showed that an 
appropriately centered and standardized uniformly distributed vector on the set A1„(M) con- 
verges weakl y to a multivariate normal di s tribut ion. This work was continu ed and substantially 
extended by iGamboa and Lozada-Chang] (120041 ) and iLozada-Changl (120051 ) , who derived a cor- 
responding large deviation principle in the one- dimensional case. 

In the present paper we will investigate related questions for the moment space (11.21) corre- 
sponding to the matrix measures on the interval [0, 1]. More precisely, we consider a uniformly 
distributed vector {Si^n, ■ ■ ■ , Sn,n)'^ on A^„(M) C (5p(M))" (for a precise definition see Section 2) 
and show that the vector of the first k matrices converges weakly after an appropriate standard- 
ization, that is 



V^{A-'^Ip){Si,n-S^„...,S, 



k,n 



Sir 



V 
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where A G R^^'^ is a matrix which will be specified in Section 2, ® denotes the Kronecker 
product, Sj = s^Ip, 



are the moments of the arcsine distribution and Gi, . . . ,Gk are independent random pxp matri- 
ces, each distributed as the Gaussian orthogonal ensemble. The proof is based on the introduction 
of new "coordinates" for the moment space A1„(M). More precisely, we define a one to one map 
from the interior of 7Vi„(R) onto the product space (Op, /p)", where Op is the pxp matrix with 
vanishing entries and the open interval with respect to the Loewner ordering is defined by 

(1.4) (Op, Jp) = {Ae Sp{R) \ Op<A<Ip} . 

That is, (Op, Jp) denotes the set of all positive definite matrices A G >Sp( R) for which /p — A 



i s pos itive definite. The new coordinates are called canonical moments [see iDette and Studden 



( 120021 )]. and they are related to the V erblunsky coeffi ci ents, which have been discussed for ma- 



trix measures on the unit circle [see iDamanik et al.l (120081 )]. We show that for a uniformly 



distributed vector on the nth moment space Ain(R) the corresponding canonical moments are 
independent and have multivariate pxp Beta distributions. Each canonical moment converges 
weakly (after centering and standardizing it appropriately) to the Gaussian orthogonal ensem- 
ble, and this result will be used to obtain a corresponding asymptotic result for the vector 

'\A^(*S'l,n Si, . . . , Sk^n S^)"^ ■ 

The remaining part of this paper is organized as follows. In Section 2 we introduce the basic 
notation, define a uniform distribution on the moment space A^„(M) and state our main result. 
We also determine the volume of A^n(M) defined by ( II. 2p . In particular, it is shown that the 
volume behaves asymptotically as 2~"'^^^^^^^^^, which means that the moment space 7V1„(M) 
defines a very small part of (5p(M))". Canonical moments of matrix measures on the interval 
[0, 1] are introduced in Section 3. The proof of our main result is given in Section 4, which contains 
several results which are of own interest. In particular we prove the weak convergence of the 
(appropriately standardized) multivariate Beta distribution to the Gaussian ensemble. Section 
5 extends these results to random moment sequences corresponding to matrix measures with 
complex entries. Roughly speaking, a corresponding weak convergence result is still available, 
where the Gaussian orthogonal ensemble has to be replaced by the Gaussian unitary ensemble. 
Finally, the proofs of some technical results are deferred to an Appendix in Section 6. 



2 The uniform distribution on the moment space of ma- 
trix measures 

Throughout this paper let (iSp(M), i3(iSp(M))) denote the measurable set of all symmetric pxp 
matrices with real entries, where jB(5p(M)) is the Borel field corresponding to the Frobenius norm 
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I |y4| I = A/tr(y42) on 5p(]R). In order to define a uniform distribution on tlie matrix moment space 
7V1„(M) we consider on 5p(M) tlie integration operator 



:2.i) 



dX = Y\ dxij , 



the product Lebesgue measure witli respect to tlie independent entries of a symmetric matrix. 
For an integrable function / : 4Sp(M) — > M the integral 



[2.2] 



f{X)dX 



is thus the iterated integra l with respect to each element Xjj, i < j [see e.g. iMuirheadI (Il982l ) 



or 



Gupta and Nagan (120001 )]. We will repeatedly integrate functions F : 



>Sp(M), in this 



case we define 
(2.3) 



F{X)dX 



{F{X)\,dX 



It was shown in iDette and StuddenI (120021 ) that the moment space A1„(M) is compact and has 
non empty interior, say Int(A^„(]R)), which enables us to define a uniform distribution on Ain{ 
To be precise we introduce the matrix valued Hankel matrices 



^2.4) H 



■2m 



and 



(2.5) H 



2m+l 



( So ■ 
( Si 

\ Sm+1 



Sm ^ 
S2m ) 

Sm+1 ^ 
S2m+1 ) 



H 



2m 



H 



^ Si — S2 ■ ■ ■ Sm ~ Sm+1 ^ 
\ Sm Sm+1 ■ ■ ■ S2m~l S2m ) 

( So ~ Si ■ ■ ■ Sm ~ Sm+1 ^ 
\ Sm Sm+1 ■ ■ ■ S2m S2m+1 / 



2m+l 



Dette and StuddenI (120021 ) showed that the point {Si, . . . , Sn)'^ is in the interior of the moment 
space Ain{^) if and only if the matrices H„ and if„ are positive definite. 
For a point {Si, . . . , Sn)'^ G 7W„(M) we define 





= {Sm+1 


" " ■ S2m) 




L±2m-1 


= {Sm: ' 


■ 5 S2m-l) 






= {Sm ~ 


Sm+1} ■ ■ ■ ; S2m-1 


~ S2m) 


"'2m-l 


= {Sm 


Sm+1} ■ ■ ■ ) S2m-2 


— S2m-l) 
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and consider the p x p matrices 

(2-6) 5^+1 = hlH-\h„, n>l, 

(2.7) = Sn-hlH-\K, n>2, 

(for the sake of completeness we also define = and = Ip, S2 = Si). Note that S~_^_'^ 
and 5*^+1 are continuous functions of the moments Si, . . . , Sn and that S^^i < S^_^_i if and only 
if (^1, . . . , Sn)^ G Int(A<„(M)). Moreover it follows that 

(2.8) Mnim = {{Si,...,Snf \ Si <Si<St,...,S- <Sn<S+}, 

(2.9) Int(A<„(M)) = {{Si,...,Snf\Si <Si<St,...,S- <Sn<S^}^(lS 



see 



Dette and StuddenI (120021 ) for more details]. Consequently the |np(p+l)-dimensional volume 



of the moment space 

(2.10) ViMnim = [ dSi...dSn, 

defined by n iterated integrals of the form as in (12. 2p is positive. A uniform distribution on 
7V1„(M) is well defined by the density 

(2-11) f{Si, . . .,Sn) = ■ ■ - ySn) ■ 

For the sake of brevity we use the notation {Si, . . . , Sn)'^ ~ W(A1„(M)) throughout this paper. 
The following result gives the volume of the moment space A^„(R). The proof will be given in 
Section 3 where more powerful tools have been developed for this purpose [see Remark 3.6]. 

Theorem 2.1. For the real moment space A^„(M) defined in (11.21) we have 

n 

(2.12) V(X„(M)) = n^p(i^(P+l)>l^(P+l)) , 

k=l 

where Bp{a, b) denotes the multivariate Beta function 

(2.13) Bp{a,b) ■.= ^^lll^^^l^^ a,h>\{p-l). 
and Tp{a) the multivariate Gamma function 



Tp{a) : 



detX"-(^'+i)/2g-*^(^)dX 

JJr(a-i(z-l)) , a>\{p-l). 



P 

P(P"l)/4 



1=1 
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As a simple consequence of Theorem 2.1 we obtain by Stirling's formula the following approxi- 
mation for the volume of the nth moment space 

logV(.M„(R)) 
lim --— = 1 , 

which shows that A1„(R) consists only of a very small part of (5p(R))"'. We will conclude this 
section with the main result of this paper, which gives the asymptotic distribution of the vector 
of the first k components of a uniform distribution on A4n(R)- For this purpose recall that a 
random symmetric matrix X is governed by the Gaussian orthogonal ensemble (GOE), if its 
density is given by 



(2.14) /(X) = (27r)-*'/V-P(P-^)/^ 



e 



Theorem 2.2. // Sn,n = {Si^n, ■ ■ ■ , Sn,n)^ ~ W(-^n(I^)); then an appropriate standardization 
of the vector Sk,n = (>S'i,n, • • • , Sk,n)'^ converges weakly to a vector of independent Gaussian 
orthogonal ensembles, that is 

v/4n(p+l)(A-i ® /p)(Sk,n - S°) G . 

Here S° = (s?Jp, . . . , s^Ip)'^ , denotes the mth moment of the arcsine distribution on the 
interval [0, 1] defined in A is a k x k lower triangular matrix with elements aij defined by 

(2.15) «M =2''^^'(^^'^-) J<'^ 

and G = (Gi, . . . , Gk)^ G (iSp(]R))^ is a vector of k independent Gaussian orthogonal ensembles, 
i.e. Gi,...,Gk i.i.d. ~ GOE. 

The proof of Theorem 2.2 is complicated and given in Section 4, which contains also several results 
of independent interest. It requires some explanation of the relation between the ordinary and 
canonical moments of a matrix measure, which will be presented in the following section. 



3 Symmetric canonical moments of matrix measures 

Let (5*1, . . . , Sn)'^ € Int(A^„(R)) be a vector of moments of a real matrix measure on the interval 
[0,1] and recall that the matrices and 5*^ given in (12. 7p and (12.61) . respectively, depend 
only on the moments 5*1, .. . Sk-i- The corresponding canonical moments of the moment point 
(5*1, . . . , Sn)"^ are defined by 

(3.1) U, = {St - S^)-\S, -S^), k = l,...,n, 
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whenever 5*^ — 5*^ > Op, otherwise they are left undefined [see lDette and StuddenI (120021 )] . Note 
that in general the matrices Ui, . . . ,Un are not symmetric and a symmetric version of canonical 
moments, say {Ui, . . . , Un)'^, can easily be obtained by the transformation 



(3.2) 



s^r'/\s,-s-)is^-s-) 



-1/2 



Throughout this paper we will work with both definitions of these matrices (no t e tha t the 
matrices Uk and Uk are similar). It is shown in Theorem 2.7 in lDette and StuddenI (120021 ) that 
the "range" of the fcth moment can be expressed in terms of the canonical moments Ui and 
Vi = I„ — Ui for i = 1, k — 1: 



(3.3) 



St-Sr = U,Vi...Uk-iV, 



k-l 



k 



, n 



Furthermore, the moment space A1„(M) is convex and by the discussion in Section 2 it follows that 
A1n(K) has non empty interior. Consequently, any random variable {Si, . . . , SnY ~ W(A1„(R)) 
with density (12. lip satisfies P{{Si, . . . , Sn)^ € Int(7V(„(]R))) = 1 and the corresponding canonical 
moments f/i, . . . , f/„ and f/i, . . . , f/„ are well defined with probability 1. Moreover, it is easy to 
see that < Sk < 5*^ implies Op < f/fc < Ip whenever (5*1, . . . , SnY ^ Int(A^„(]R)), since the 
Loewner ordering is not changed by pre and post multiplication with positive definite matrices. 
Given 5*1, .. . 5*^-1, the moment 5*^ can be calculated from the canonical moment U^ and therefore 
equation (13. 2p defines a one to one mapping 



(3.4) 



Int(;W„(M)) 

(5*1, . . . , S-nf^ 



(Pp{Si 



.Sn) = {Uu...,Unf 



from the interior of the moment space onto the "cube" (Op, Jp)" defined by (11.41) . In the following 
Lemma we collect some interesting properties of the matrix valued canonical moments, which 
will be useful in the following discussion. The proof can be found in the Appendix. 



Lemma 3.1. 



(a) If fi is a matrix measure on the interval [0,1] with corresponding canonical moments U^ 
and u = fi"' is the measure induced on the interval [a, b] by the transformation 7(x) = 
{b — a)x + a {a < b) with corresponding canonical moments f/^, then 

U" = U^ 

whenever the canonical moments are defined. In other words: the canonical moments are 
invariant under linear transformations. 

(b ) If the matrix measure is symmetric, then 

U2n-l = \lp) 

whenever the canonical moments are defined. 
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(c) Let fi denote a matrix measure on the interval [0, 1] with canonical moments and define 
a as the symmetric matrix measure on the interval [—1,1] induced by the transformation 

(3.5) a{[-x,x])=fi{[0,x']) 
with corresponding canonical moments U^, then 

(3.6) f/^Vi = lip , u^n = ua, 

whenever the canonical moments are defined. 

The following result shows that the ordinary moments of a matrix measure can be calculated 
recursively frorn the canonical moments Uj. A similar result in the scalar case was shown by 



Skibinskvl ([19681) 



Theorem 3.2. For a moment point {Si, . . . , Sn)'^ € Int(A^ri(M)) with corresponding canonical 
moments Ui, . . . , Un, define Co = Op, Ci = Ui and 

(3.7) Q = V,^iUj, j = 2,...,n. 
Then we have 

Sn Gfi,ni 

where {Gij, i,j G (1, . . . ,n)} denotes an array of p x p matrices defined by Gtj = Op if i > j, 
Go,j = Ip and recursively by 

(3.8) Gjj = Gjj-i + Cj-i+iG'i-ij , 
whenever j > i > 1. In particular, we have 

^ ^ Cin ■ • • • ■ Ol ! 

where the index set is defined by I = {(zi, . . . ,in)\ik £ {1, • • • , n}, ii = 1, ik < ik-i + 1} ■ 
Proof: We consider the (infinite dimensional) block Hankel matrix 

(3.9) M = (5i+,).,,>o , 

which contains the moments of the matrix measure /x. Let {Pn{.x)}n>o denote the sequence of 
monic (this means that Pn{x) has leading term x'^Ip) orthogonal matrix polynomials with respect 
to /i, that is 



(3.10) j P^{x)dii{x)Pl{x) = I 



_ / Op G \in^m 
Dn G iin = m 



It was shown by lSinap and van Asschd (119941 ) that these polynomials satisfy a three term recur- 
rence relation 



Poix) 



Pi (x) = xip — Ai 



xPn{x) = Pn+l{x) + An+lPn{x) + Bn+lPn-l{x) 

We define P(x) = {P^{x), Pf (x), P^{x),...Y, and 



n > 1 



(3.11) 



Ml ip Op ■■■ 

B2 A2 Ip Op 

0„ B-i A-i In 



\ 



then the recursion can be rewritten in the form 

(3.12) JP(x) = xP{x) . 

If F{x) = {Ip, xIp, x'^Ip, ■ ■ denotes the vector of matrix valued monomials, then it follows that 

(3.13) P(x) = LF(x) , 

where L is a lower triangular block matrix containing the coefficients of the matrix polynomials 
Pn{x). The following Lemmata are proved in the Appendix. The first specifies the inverse of the 
matrix L. 

Lemma 3.3. The matrix L defined in (13.131) is non singular and its inverse K := L^-*^ is defined 
by 



(3.14) 

where the matrix R is given by 



R 



RK = KJ , 



^Op Ip Op 

Op Ip Op 



\ 



J 



Moreover K a lower triangular block matrix with the matrices Ip on the diagonal. // D = 
diag{Do, Di, . . .) is the block diagonal matrix with entries Dj defined by (I3.10p . then the Hankel 
matrix M defined by (13. 9p has the representation 



(3.15) 



M = KDK^ 
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Lemma 3.4. // fi denotes a matrix measure on the interval [0, 1] and a the corresponding 
symmetric measure on the interval [—1, 1] defined by (13.51) . then the monic orthogonal matrix 
polynomials {Pn{x)}n>o with respect to the matrix measure a satisfy the recurrence relations 

Po{x) = Ip, Pi{x) = xip, 

(3.16) xPn{x) = Pn+l{x) + CPn-l{x), n>l , 

where (n = Vn~iUn o,nd Un are the canonical moments of the measure /i. 

By Lemma 3.4 the orthogonal polynomials Pn{x) with respect to the matrix measure a on the 
interval [—1, 1] are even (if n is even) or odd (if n is odd) functions. Consequently it follows from 
Lemma 3.3 and the representation (13.131) for the block Kij G R^^^ in the position (i, j) of the 
matrix K = L~^ that Ki^j = Op if i + j is odd. Moreover, for the elements of the corresponding 
matrix J in (13.111) we have A^+i = Op and -Bn+i = Cn ? where Cn corresponds to the matrix 
measure fi. Observing (13.141) we obtain the recursion 

(3.17) Ki^2j,i = Ki^2j-l,i~l + Ki^2j-l,i+lCj+l ■ 

With the definiton 

{a m > n we define = Op) one easily sees that the matrices Gm,n satisfy the recursion (13.81) . 
Finally the representation (13.151) yields for the moments and of the matrix measures /i 
and a the relation 



— ^2n — ^2n,0 — -f^2n,0-Do-ft'o,0 — -f^2n,0 — Gn 



,n } 



where Mij denotes the p x p matrix in the position of the matrix M corresponding to 
the matrix measure a. This proves the first part of Theorem 3.2. The remaining statement is 
obvious. □ 

In the following we will study the distribution of the canonical moments corresponding to a 
random moment vector uniformly distributed on Aini 



Theorem 3.5. // (5*1, . . . , Sn)'^ ~ U{Ain{^)) is a random vector with a uniform distribution 
on the nth moment space A^„(R), then the distribution of the corresponding random vector 
of matrix valued canonical moments {Ui, . . . ,Un)'^ is absolute continuous with respect to the 
Lebesgue measure and its density is given by 

1 " 

(3.18) /(f/i,...,t/n) = viMnim Ui ' 
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where V(A^„(M)) is defined m fIZTU]) . 



Note that Theorem 3.5 shows that the random variables [/i,...,f/„ corresponding to a ran- 
dom moment vector (5*1, . . . , SnY ~ W(A^„(M)) are independent and have a multivariate Beta 
distribution, that is 

Uk ~ Betap{^{n-k + l){j)+l),\{n-k + l){j)+l)) , 

where the density of a random variable X ~ Betap{a, b) with a matrix valued Beta distribution 
with parameters a and b is given by 



(3.19) 



/(X)=5,(a,6)-^(detX)'^-(^+i)/2(det(/-X))''-(^+^)/2^(o„,,)(X) 



see 



Olkin and RubinI (119641 ) or lMuirheadI (119821 )] and the normalizing constant Bp{a, b) is defined 



in (I2.13p . Note that this definition requires a,6>(p — l)/2. 

Proof of Theorem 3.5: We first calculate the Jacobi determinant of the mapping cpp in (13. 4p 
from the ordinary to the canonical moments, which we denote by J{fp)- By the transformation 
formula (13.21) we can write (pp{Si, . . . , Sn) = {ip^\Si), . . . , (p^\Sn)), where 

: {S^,S^)-^{Op,Ip) , 

^f\Sk) = [St - s,)-'/'s,{st - 5, )-v2 _ (5+ - s,)-'/'s,{st - s,)-'/' . 

(k) 

Note that the transformation ipp is one to one and depends only on the moments Si, ... , Sk-i. 
This implies that the Jacobian J{(pp) is the product of the Jacobians of the transformations (p^^ . 
For fixed nonsingular matrices A and B with B G 5p(]R) the Jacobian of the transformation 
X I— > AXA^ + S is equal to (det A)^"*"^, see Theorem 2.1.6 in iMuirheadI (119821 ). Consequently 
we obtain with the aid of equality (13.31) 



k=l k=l 
n 

l[det{UiVi...Uk-iVk 

k=2 

n ^^— 1 



'k~i 



^(P+l)/2 



k=2 



k=l 



This gives for the density of the vector {Ui, . . . UnY 

1 



n-l 



v{Mnm) 



J{(?7i, . . . , Unf e ^p{lni{Mnm))] n det(f/fcV^fc)(" 



-fe)(p+l)/2 



k=l 
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where v?p(Int(A^„(M))) = (0^,/^)-. 



□ 



Remark 3.6. Note that the proof of Theorem 3.5 provides also a proof of the formula for the 
volume of the nth moment space in Theorem 2.1, because 

„ n— 1 

ViMnim = / lldet{Uk{Ip-UkY''-'^^'^+'^/^dUi,...dUn 

n 

= J]5p(iA;(p+l),ifc(p+l)) . 

k=l 



4 The multivariate Beta distribution and a proof of The- 
orem 2.2 

The proof of Theorem 2.2 is separated in two steps. First we investigate the asymptotic properties 
of the multivariate Beta distribution (Section 4.1). In particular, we show that a standardized 
version of the random matrix X„ ~ Betap{an, an) converges in distribution to the GOE if the 
parameter a„ tends to infinity. From this result and Theorem 3.5 we obtain a weak convergence of 
the vector Uk,n = (Ui^n, ■ ■ ■ , Uk^n)'^ of canonical moments corresponding to the first k components 
of a vector Sn,n = (•S'l,™, . . . , Sn,nY ~ if n tends to infinity. 

Secondly, we use the relation between ordinary and canonical moments of matrix measures on 
the interval [0, 1] to prove a corresponding statement regarding the weak convergence of the 
vector Sk,n = (5'i,„, . . . , Sk,nY (Section 4.2). 



4.1 Some properties of the multivariate Beta distribution 

By Theorem 3.5 the multivariate Beta distribution will play a particular role in the analysis of 
random moment sequences of matrix measures on the interval [0, 1]. This distribution on iSp(M) 
can easily be defined by its density (13.191) . Since the density depends on X only through the 
determinant of X or Jp — X, the distribution of a multivariate Beta distributed random variable 
X is invariant under the transformation X 0X0^ for any orthogonal matrix O G 0{p), 
where 

0{p) = {0e M^^^l 00^ = Jp} 



denotes the orthogonal group. For some properties following from this invariance see I Gupta and Nagar 



( I2OOOI ). chapter 9.5. The eigenvalues of a multivariate Beta distributed random variable follow 



the law of the Jacobi ensemble. To be precise recall that the Jacobi ensemble is defined as the 
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distribution of a vector A = (Ai, . . . , Xp)^ with density 



(4.1) 



(0,1) 



(A. 



i=l 



where A(A) = ni<j(Ai ~ Xi) is the Vandermonde determinant, a,b, (3 > and the constant cj is 
given by 



(4.2) 



A r(i + f)r(a + & + f(j9 + j-2)) 

l\ Til + 4 



+ fj)r(« + f(j-i))r(6+f(j-i)) 



see for example iDumitriu and EdelmanI (120021 ) . For the sake of simphcity we write 



(4.3) 



Xr^J, 



(a,fe) 



/3 



if a random vector A = (Ai, . . . , Xp)'^ has density (14.11) . Usually only the cases P = 1,2 and 4 
are conside r ed cor responding to matrices with real, complex and quaternion entries, respectively 
[see iDysonI (119621 )]. A symmetric random variable X ~ Betap{a,b) can be factorized as X = 
Odiag{X)0^ , where O G 0{p) and diag{X) is a diagonal matrix containing the eigenvalues 
of X. Integration with respect to the orthogonal matrix O shows that the eigenvalues are 
distribute d acco rding to the Jacobi ensemble j7'^('^~(p~^)/^'^~(p~^)/^)^ for the calculation we refer to 
MuirheadI (119821 ). We now make use of the invariance of the multivariate Beta distribution and 
the distribution of the eigenvalues and calculate the first moments 

(4.4) E[X^] = j X^f{X)dX 

of a multivariate Beta distribution. 

Lemma 4.1. Suppose X ~ Betap{a,b) , then the moments defined by (14.41) satisfy 



(4.5) 



E[X''] = cjp , 



where the constant G M depends on the parameters a, b and p. In particular, we obtain for the 
first two moments of X 



(4.6) 
(4.7) 



E[X] 
E[X2] 



a + b 



(a + 6)(a + 6+ 1) 



a + 1 + (p - 1) 



2a + 26 - 1 ' 



Proof: Because of the invariance of the multivariate Beta distribution we obtain for any orthog- 
onal matrix U 

(4.8) E[X'']U = E[UX''U'^]U = UE[X''] . 
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Therefore E[X'^] commutes with all orthogonal matrices, which gives E[X'^'] = c^/p. The real 
constant Ck can be determined by 



(4.9) 



1 



1. 



Ck = -trE[X'=] = -E[trX'^] = -E[A^ + ■ ■ ■ + A, 



p p p 

where the distribution of the eigenvalues Ai, . . . , Ap is the Jacobi ensemble with paramters a — 
|(p — 1), 6 — |(p — 1) and P = 1. Therefore the moment E[X'^] is given by 

(4.10) E[At] ■Jp = c, [\.. /"'At|A(A)inA;-(^+^)/'(l-A,)^-(^+i)/2rfA - I,. 

Jo Jo 



This integral is known as Aomoto's generalization of the Selberg-integral [see lAomotd (119881 )]. 

'2-y 



Aomoto showed that the eigenvalues of the Jacobi ensemble satisfy 



(4.11) 



E[Ai • . . . ■ A, 



n 



a + 7(p - i) 



\a + (3 + 7(2p 



for 1 < m < p. By a similar method iMehtal (120041 ) gets the recursion 
(4.12) {a + (3 + l + 27(p - l))E[Xl] = {a + 1 + 2^{p - l))E[Ai] - -f{p - l)E[AiA2] 

We combine equation (14.111) and (I4.12p and obtain 

a + 7(p - 1) 



(4.13) 
(4.14) 



E[Xi] 



{a + P + 27(j9 - l))(a + /? + 1 + 27(^9 - 1)) 
X f(a + l+7(p-l))+7(p-l) /5 + 7(p-l) 



a + P + 2-f{p - 1) - 7 



This completes the proof of Lemma 4.1 if we set a = a — ^(p — 1), P = b — ^{p — 1) and 7 = |- 



□ 



As a next step we state a result concerning the asymptotic properties of the multivariate Beta 
distribution if the parameters tend to infinity. 

Theorem 4.2. Assume that ~ Betap{an, an) for a sequence a-n/n 'j E M+, then 

(i) Xn ^ ' 

n— >oo 

(it) V8^(X„-i/p) 

n^oo 

where the random variable G is distributed according to the GOE. 
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Proof: Let || ■ || denote the Frobenius norm, then we obtain by Lemma 4.1 
E|||.f„ - = trE[(A'„ - E[X,,]f] = tr(EKl " E[X„]2 



a„ + 4 \ 4a„ — 1 / n^oo 



The proof of (ii) is based on the convergence theorem of IScheffe (119471 ). by which it suffices 
to show that the density fn of the standardized random variable y/8'jn{Xn — |/p) converges 
pointwise to the density / of the GOE given in fl2.14p . The density is given by 



MX) =i?;i(a„,a„)(87n)-P(^+^)/^ 



a„-(p+l)/2 ^ N a„-(p+l)/2 

X det ( + ^/p ) det ( -Jp - ^==X Ii-V2^i„V2^i,){X) 



/ \ a„-(p+l)/2 

We can diagonahze a fixed matrix X G iSp(]R) as X = Odiag{X)0^ , where O G 0{p) is an 
orthogonal matrix and A = (Ai, . . . Xp)^ are the eigenvalues of X. Therefore is easy to see that 
each factor in the last formula satisfies 

^ a„-(p+l)/2 

-X- 

2 

an-{p+l)/2 



P / Na„-(p+l)/2 



n— >oo 

i=l 



As n tends to infinity, we obtain by Stirling's formula 



,(8^^)-p(p+l)/42-2pa„+p(p+l)£p(^^ 



^^-p(p-l)/4(^g^^^-p(p+l)/42-2pa„+p(p+l) J-j- 



P r(2a„-i(z-l)) 



r(a„- 1(2-1)) 



2 



^^-p(p-i)/4^27r)-P/2(i + o(i)) 



In other words, the normalization constant of the density /„ converges to the normalization 
constant of the GOE, which completes the proof. □ 
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By Theorem 3.5 it follows that for a vector of matrix-valued moments Sn,n = . . . , S'„^„)^ ~ 
h({A4n{^)) chosen uniformly from the moment space Ain{^) the corresponding canonical mo- 
ments Ui^n, ■ ■ ■ , Un,n are independent multivariate Beta distributed. As n tends to infinity the 
parameters of the Beta distributions behave as The following Theorem is thus a direct 

consequence of Theorem 4.2 with 7 = |(p -|- 1). 

Theorem 4.3. Assume that Sn,n = ('S'i,„, . . . , Sn,nY ~ W(A^n(M)) and let Uk,n = (f^i,n, • • • , Uk,n)'^ 
denote the vector of the first k canonical moments corresponding to the random variable Sn,n- 
Then 

V4(p+l)r2 (Uk,n - U°) Gk , 

n^oo 

where U° = |(/p, . . . , Ip)'^ and Gk consists of k independent matrices of the GOE. 

It follows from calculations in the scalar case that the canonical moments of the arcsine distri- 



bution defined in (11. 3p are all equal 1/2 [Skibinskyl (119691 )]. Therefore we obtain 



U° = -(I I 



T 

2 V-'P' ■ ■ ■ j-^pj ■ 

^p{{{s°Jp, . . . , slip) 



= MK) 

In other words the vector U° used in the centering of Theorem 4.3 contains the canonical 
moments corresponding to the matrix measure /i defined by 

(4.15) dfi{x) = — -.Ipdx . 

'n-y'x{l — x) 

For this reason the sequence of moments S° of the matrix measure defined by (14.151) can be 
viewed as the "center" of the moment space Aij 



4.2 Asymptotic properties of random moments 

In this Section we will use the results of Section 4.1 to prove Theorem 2.2. The basic idea of the 
proof consists of two steps. First we show that the inverse of the mapping 

(4.16) ipp : Int(A^,(M)) — ^ (0^, Ip)" 

defined as in (13.41) is different iable in a sense defined below, secondly we use this property and 
Theorem 4.3 to establish the weak convergence of the vector Sk,n of the first k components of 
Sn,n = (5'i,n, • • • , Sn,n)'^ ~ W(7V1„(M)). For this purposc recall that 

(4.17) v^(Sk,„ - s°) = Mv;\^Kn) - ^;'(u°)) . 
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where Uk,n = (f^i.n, • • • , Uk^n)'^ denotes the vector of canonical moments corresponding to Sk,n 
and U° = |(/p, . . . ,IpY ■ In the scalar case p = 1 the quantity in ( 14.17P can be reduced by 
differentiating the mapping that is 

(4.18) v^(sk,„ - s°) = v^((/^r'(uk,n) - Vl\0) 

d 

= V^^-^(Uk)(uk.n -u^) + y/n o(||uk,n - u°) , 

where for the sake of readability the lower capital symbols Sk,n, s^, Uk,n and u° denote the 
moment vectors Sk.n, S^, Uk,n and in the case p = 1. Note that 



u") = A 



where the elements of the matrix A were found by lChang et al.l (119931 ) and are defined by (12.151) . 
In order to study the general matrix case we introduce the following concept of differentiability. 

Definition 4.4. Assume that S C (5p(M))" is an open set. A mapping ^ : S (]^p>^p^"i 
called matrix differentiable in a point M° G S, if there exists a matrix L G such that 

(4.19) $(M° + H) - $(M°) = LH + o(| |H| |) . 

In this case the matrix derivative of ^ at the point M° is defined by ^(M°) := L . 

Note that matrix differentiability is a stronger concept than total differentiability and that a 
linear mapping $i(M) = AM + B is matrix differentiable with |^ = A. On the other hand the 
mapping $2(^) = is only matrix differentiable at the points = mlp. It is easy to see 
that matrix differentiability has the usual properties and we note for later reference 



(4.20) 



9$ 9$ 



dM \dM ' ' dM 
if $ = ($^, . . . , $^)"^ is matrix differentiable, and 

(4.21) V M°) = ^(M°)^(M°) + $(M°)-^(M°) . 

if m = 1 , $ and are matrix differentiable in M° and \E'(M°) = c/p with c G M. Our next 
result shows that the inverse of the mapping yjp defined in (13. 4p is matrix differentiable and gives 
the derivative. The proof is complicated and given at the end of this Section. 

Theorem 4.5. The mapping (p~^ : {Op,Ip)^ lnt{Aik{M.)) defined by (13.41) is matrix differen- 
tiable at the point U° = |(/p, . . . , Ip)^ with 

(4.22) ^(U°) = A®/p, 
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where A is the lower triangular matrix defined in fl2.15p . 



With the aid of Theorem 4.5 we are now in a position to complete the proof of Theorem 2.2. 
More precisely we obtain from (14.171) and fl4.22p 

V^(Sk,„ - S°) =v^(¥., ^(Uk,n) - 

= V^{A ® /p)(Uk,n - U°) + Op(| |Uk,n - U°| I) 

and the assertion of Theorem 2.2 follows because (14.61) and (14. 7p yield that the expectation of 
'^llUk.n — U°|p converges to pk/8, which implies that n||Uk,n — U^|p = Op(l). Also note that 
A® Ip is nonsingular and {A Ip)~^ = ® Ip- □ 



Proof of Theorem 4.5: We first study for 1 < m < k the mapping 



(4.23) V 



{Ui,...,Uk)^U„ 



where U = {Ui, . . . , Uk)^ G (Op, Ip)'' is a vector of (symmetric) canonical moments defined by 
(13. 2p and Um the mth non symmetric canonical moment defined by (13. ip . Note that Um = 
Dm^'^UmDll'^ where 

Dm = -Dm(U) = — , 



and D m satisfies the recursion Dm+i = Dli'^UmVmDm^ , Di = Ip [see Theorem 2.7 in Dette and Studden 
(I2OO2I )]. Obviously D^ depends continuously on Ui,...,Um~i- At the point U° we have 
D„(U«) = An(|/p, • • • , = and V(U°) = lip. With the notation = /^™(U° + 

H) we obtain for H = (Hi, . . . , Hk)^ e 5p(M)^ 



V^(U° + H) - ^(U°) = ^„^/^(i/p + i/„.)^;t/' - i/p 

= ^v^ra + {D^^'^IIraD\l^ " -f^»7 



The remainder can be estimated as follows 

= b~^/'^H i b^'^ - D^'^) + i b-^/'^ -DH 

= o(l|H||). 

This yields l^(U'') = e^® Ip and as a consequence 



^(U°) = ^(u°)®/. 



p 1 
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where as in fl4.18p Um denotes the mth component of the vector u of canonical moments in 
the case p = 1 and u° is the vector of scalar canonical moments corresponding to the arcsine 
distribution in (11.31) . A similar argument shows that 

where Vm = ^ — ^m- By (14.211) products of canonical moments Um and Vm are also matrix 
differentiable at the point U° (for sums this statement is trivial) and the derivative is the 
Kronecker product of the derivative in the case p = 1 with the unit matrix. For example 
we can calculate for m ^ 1 

^Cm/j-rON _ 9Vm-lUm /ttO\ _ It 9Vm-l /ttO\ , 1 r ^J^(JjO\ 
\^ ' nxT \^ } — 2^P nxT ^ ' ^ 2^P nxT 



[em-em-i) ®Ip = 7. (u)0/p 

ou. 



and for m 



p 



Finally Theorem 3.2 shows that the mth moment Sm is equal to the sum over products of the 
matrices Cm- Therefore the matrix derivative of Sm with respect the canonical moments U is 
given by ^(U") = ^(u°) ® Jp and fOOD yields 



^(U°) = ^(u°)®/, = ^<S/,. 
This completes the proof (note that A® Ipis non singular). □ 



5 Complex random moments 



To a large extend, the case of complex matrix measures can be treated analogously to the case 
of real matrix measures. For the sake of brevity we only state the results in this Section and 
omit the proofs. The fcth moment of a complex matrix measure on the interval [0, 1] is defined 

as 

(5.1) Sk= ! x^d^{x) eSp{C)] A; = 0,1,2,... 

Jo 

where iSp(C) denotes the space of p x p hermitian matrices. The complex nth moment space 



(5.2) 



A^„(C)= (5i,...,5, 



S, 



xM/i(x), J = 1, . . . , n ^ C (>Sp(C)) 
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is characterised by the equations fl2.8l) as well [see Dette and StuddenI ( 2002 )]. Here A* = 
denotes the conjugate transpose of the matrix A. For a point {Si, . . . , Sn)* G Int(A1n(C)) the 
complex canonical moments Ui, . . . ,Un are therefore well-defined, where as in the real case 



(5.3) 



and the hermitian matrices 5^ and are defined as in fl2.6p and (12.71) . respectively. The 
integration operator changes on Sp{C) to 



(5.4) 



dX = dxii dReXijdlmx, 



i=l 



i<j 



that is, we integrate with respect to the independent real entries of a hermitian matrix. 
Note that in this case for a nonsingular matrix A G Sp{C) the Jacobian of the transformation 
X AX A is given by (det A)^^. The law of a random variable X G Sp{C) is called complex 
multivariate Beta distribution with parameters a, 6 > p — 1 if its density is given by 



(5.5) 



/(X) = (5f (a,6))-MetX»-Pdet(/,-X)''-^/(o„,,)(X) 



('2) 

and we denote this property by X ~ Betap {a,b). The normalizing constant is the complex 
multivariate Beta function 



r?\a)r?\b) 

(a + 6) 



where T^\a) 



p(p-l)/2 



n f-i r(a — i + 1 ). F or a more genera l discu ssion of the complex 
Beta distribution we refer to iKhatril (119651 ) and iPillai and Jourid (jl97ll ). The eigenvalues of 



(2) 

a B etttp (a, &)-distrib u ted r andom variable follow the law of the Jacobi ensemble J'r 



{a-p+l,b~p+l) 



see 



Pillai and Jourid (Il97ll )] and similar arguments as given in the proof of Lemma 4.1 show 
that the first moments of a random variable X ~ Beta^\a, b) are given by 



E[X] 



a + b 



a + 1 + (p - 1) 



{a + b){a + b+l) \" ' ~ ' ~'a + b-l^ 
Proceeding as in Section 3, we get the following result for complex canonical moments. 

Theorem 5.1. Let Sn.n = {Si^ni ■ ■ ■ ■, Sn.n)* be uniformly distributed on the complex moment 
space Aini'C) defined in (15. 2p . then the corresponding canonical moments Ui^n, ■ ■ ■ , Un,n in- 
dependent and for k = 1, . . . ,n Uk^n is complex multivariate Beta distributed with parameters 
{p{n — k + l),p{n — k + 1)). 

('2) 

For a sequence of complex random variables X„ ~ Betap {an, an) an analoge of Theorem 4.3 
holds, where in the limit the Gaussian orthogonal ensemble has to be replaced by the Gaussian 
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unitary ensemble (GUE). Recall that a, p x p hermitian matrix of the GUE is characterized by 
the density 

(5.6) f{X) = (27r)-f/V-f(^'-^)/2g-itrX2 _ 

(2,) 

Theorem 5.2. Assume that Xn ~ Betap {an, an) for a sequence an/n — > 7 G then 
(i) Xn — ^ \lp , 



(ii) y/W^{Xn-\lp) -^G , 

n^oo 

where the random variable G is distributed according to the GUE. 

The remaining arguments in Section 4 remain essentially unchanged, which yields the following 
result on the weak convergence of random complex moments. 

Theorem 5.3. // Sn,n = ('S'i,„, . . . , S'„_„)* ~ f/(A^„(C)), then the standardized vector of the 
first k moments Sk,n = {Si^n, ■ ■ ■ , S^^n)* converges weakly to a vector of independent Gaussian 
unitary ensembles, that is 



V^(A-^ /p)(Sk,n - S«) G . 

n— »oo 

The matrix A and the vector are defined as in Theorem 2.2 and G = {Gi, . . . ,Gk)*, with 
Gi,...,Gk i.i.d. ~ GUE. 



6 Appendix: Proof of auxiliary results 
6.1 Proof of Lemma 3.1 

(a) We denote by Sn and T„ the nth moment of the matrix measure ji and v, respectively, then 
a straightforward calculation yields 

n-l , X 

(6.1) ^" = E r «""(^ - «)''^^ + - «)"'^- ■ 

Note that T^ {T~) is the unique maximal (minimal) matrix with respect to the Loewner ordering 
such that for fixed To, ... , T„_i the vector (Tq, . . . , T„) is an element of the moment space of 
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the matrix measure on the interval [a, b] . Therefore we obtain (note that the specification of 
To, ... , T„_i determines Sq, . . . , Sn~i) that 



i=0 



This yields Tn — = {h — a)"'{Sn — S^), and the assertion (a) of Lemma 3.1 follows from the 
definition of the canonical moments in (13.21) . 

(b) We consider the transformation 0(x) = 1 — x and the measure v = jJ' = n. The same 
arguments as in part (a) show 

+ + + - 

which implies for the corresponding canonical moments 

(6.2) U^^ = Ul , f/^„_i = /, - . 

Because /i = i/ we obtain ?7^„„i = f/^„_i = Ip — Ul^^^^^, which yields f^2n-i — 

(c) We obtain for the moments S*^ of the matrix measure a 

(6.3) 

where Si, S2, ■ ■ ■ denote the moments of /i. The measure a is obviously symmetric and (b) yields 
^2n-i — \-^p- From (16. 3p we have for the even moments 

which yields S*^" = S~, 5*2^ = S*;]". Consequently it follows 

□ 



6.2 Proof of Lemma 3.3 and 3.4 

Proof of Lemma 3.3: From (13.131) and (I3.12p we obtain (observing that the matrix R acts as 
a shift operator) LRF(x) = xP{x) = JP(x) = JLF(x), which yields 

(6.5) LR = JL . 
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It is easy to see that the matrix L is non singular and that the inverse matrix K := L ^ is again 
a lower triangular block matrix with matrices Ip on the diagonal. From ( 16.5P we therefore obtain 

RK = KJ . 

On the other hand F(x) = KP(a;) and by the orthogonality relation (13.101) it follows 

M ^ /f(.,..(.)F-,.) ^ K . /p,.).M.)P-(.) ■ ^ KDK- . 
where the matrix D = diag{DQ, Di, . . .) is defined by fl3.10p . □ 



Proof of Lemma 3.4: It follows from Favard's theorem [see ISinap and van Asschd (119941 ) or 
Dette and StuddenI (120021 )] that there exist matrices An, Bn such that the polynomials {-Pn(a^)}n>o 



orthogonal with respect to the matrix measure a satisfy a three term recurrence relation 
Po{x) = Ip, Pi{x) = xip - Ai, 

xPn{x) = P„+l(x) + An+lPn{x) + Bn+lPn-l{x), n>l. 

We define y = |(a; + 1) and obtain from Dette and StuddenI (2002) for the monic orthogonal 

polynomials Rn{y) = '^~^Pni'^y — 1) with respect to the measure a = a 2*^^+^^ on the interval 
[0, 1] the recursion 

Ro{y) = h 



Ri{y) = ylp-C[ , 



(6.6) 



yRniy) = Rn+M + (Cl„+1 + CLf 'Rniy) + iCL-lClnVRn-liy), U > 1 



where = V^^iU^ for n > 2, = Ui and denote the canonical moments of the measure a. 
Observing Lemma 3.1 (a) and (c) it follows that Q = = ^Ip and for ri > 1 



(6.7) C2„ = V,l_,U^n = V^n-iUL 

(6.8) Cl„+i = |(/p - ^n) = |K . 
This yields for the polynomials {Pn{x)}n>o 

Po{x) = Ip, Pi{x) = xIp 

and (16.60 simplifies to 

i(x + l)2-"P„(a;) = 2-""^ (Pn+i(x) + Pn{x) + CP^ 
which proves the assertion of Lemma 3.4. 



2 ) 



-l(^)) , 



n>l 



□ 
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